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ABSTRACT 


Linear optimal regulator theory is applied to a nonlinear simulation 
of a transport aircraft performing a helical landing approach. A closed 
form expression for the quasi-steady nominal flight path is presented 
along with the method for determining the corresponding constant nominal 
control inputs. The Jacobian matrices and the weighting matrices in the 
cost functional are time varying. A method of solving for the optimal 
feedback gains is reviewed. The control system is tested on several 
alternative landing approaches using both three and six degree flight 
path angles. On each landing approach, the aircraft was subjected to 
large random initial state errors and to randomly directed crosswinds. 

The system was also tested for sensitivity to changes in the parameters 
of the aircraft and of the atmosphere. Performance of the optimal 
controller on all the three degree approaches was very good, and the 
control system proved to be reasonably insensitive to parametric uncer- 
tainties. It did not perform as well on the six degree approaches, and 
a modification to these flight paths is proposed for the purpose of 
improving performance. 
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INTRODUCTION 


Aircraft traffic in the neighborhood of commercial airports 
has been the subject of intensive study during the past several 
years. Noise, collision avoidance, airspace congestion, air 
pollution, and dangerous wing tip vortices are all problems which 
are being studied by engineers today. In the present scheme, air- 
craft are typically sequenced to follow one another down s three 
degree flight path on a straight-in approach to the runway. This 
path very often results in low altitude flying at relatively high 
power settings over residential or densely populated business 
districts which is objectionable, at least from a noise standpoint 
if not from the standpoint of safety. Recently, steeper, six 
degree flight paths and two segment approaches (Ref. 12) have been 
proposed to eliminate these problems and are currently being tested. 
This paper proposes an alternative to these landing approaches. It 
consists of descent along a helical path to a low altitude where 
the straight-in approach is intercepted and followed. This flight 
path has the advantage of keeping the aircraft at low power and at 
relatively high altitude except in the immediate area of the airport. 
The helical path is quasi-steady in that the flight path angle and 
velocity are maintained constant during descent and the control 
inputs are fixed except to null flight path errors . 
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Implementation of the proposed automatic guidance system 
requires two factors which are not at present universally available. 
The two factors are an on-board flight control computer and a 
fairly precise knowledge of the aircraft's position in space. 

However, with the advent of automatic landing systems, flight 
computers are becoming much more common on commercial transports 
and the current installation of the microwave landing system (MLS) 
equipment at airports throughout the country will make the necessary 
position data available (Ref. 12). 

The method proposed in this paper uses the theory of optimum 
control which is treated in references 1 through 5. The theory is 
applied to a three degree of freedom nonlinear aircraft simulation. 

In order to pose the problem as a linear regulator, the equations 
were linearized about a nominal state (helical path) and control 
trajectory. The resulting near optimum feedback control gains were 
then tested in the nonlinear model. This work is an extension of 
that presented in references 6 and 8. The extension includes 
increasing the number of states to six and the number of controls to 
three. The additional states removes the restriction that altitude 
be monotonically decreasing and the inclusion of thrust, autothrottle, 
as a control, allows treatment of aircraft other than gliders. 
Reference 6 concerns itself with a HL-10 lifting body reentry vehicle 
ruid reference 8 with an unpowered Boeing 707 with no external wind 
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disturbances. Another extension of this work is inclusion of the 
quasi-steady helical path with its consequent piecewise constant 
control inputs. This factor greatly reduces the problem of computer 
storage of the state and control trajectories. The method for 
computing these quasi-steady paths along viv* the method for 
computing appropriate control inputs are presented here- In addition, 
a derivation of the equations of motion is given, along with the 
method of solving for the near optimum feedback gains. It should be 
noted thpt this study, as do references 6 and 8, uses angle of attack 
and bank angle as control inputs. This is a simplification of the 
dynamic model which does not include control surfaces such as 
elevators, ailerons, or rudder. 

Feedback gains were obtained for several candidate approach 
paths including both three and six degree flight path angles and turns 
of 180, 2?0, and 360 degrees. An overhead approach was the principal 
path studied, i.e. one in which the aircraft crossed over the runway 
at a ninety degree angle heading east at an altitude of 550 meters. 

At about 1800 meters east of the runway center line it intercepted 
the helix and descended down it, making a total turn, to the right, 
cf 270 degrees. The turn war completed at an altitude of 100 meters 
where the aircraft intercepted the stkaight-in approach. At this 
time, it was about 1900 meters from the touchdown point on the runway. 
The plane continued to descend until, at twelve meter altitude, it 
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executed a flare Maneuver. The effectiveness of the control system 
vas tested in the nonlinear simulation by introducing initial errors 
in the states, cross winds , and uncertainties in the aerodynamics 
of the aircraft. 

The aircraft treated in the study is a Boeing 737-100, which is 
a two engined snail commercial transport used for relatively shcrt 
h aul flights. Langley Research Center has recently procured one of 
these aircraft from Boeing and it is being used to study and test 
advanced a utom a t ic control and avionics concepts in the neighborhood 
of the air terminal . Whereas the results presented in this paper 
shew feasibility, it is hoped that they will be applied first to a 
more sophisticated six degree of freedom simulation of the aircraft, 
and then eventually to be tested on the aircraft itself. 



A1ALYSIS 


Introduction 

The purpose of this section is to present the mathematics 
involved in simulating the aircraft and in implementing the feedback 
control law. Hot ice that in figure 1, the aircraft is modeled by a 
set of six nonlinear ordinary differential equations represented by 
the single vector equation 

^*f(x,u,t) ( 1 ) 

This set of equations was integrated nwerically using a fixed step 

fourth order Runge Kutta algorithm. The feedback control law used 

in this study requires that the nominal state trajectories, x q , 

the nominal control time histories, u , and the feedback gain 

n 

matrix, K, be stored in an on-board computer. However, it will be 
shown that by choosing particular nominal trajectories, the vector 



analytic vector function of time. This eliminates some of the 
on-board storage requirements and makes the control system simpler 
and more flexible. 

As simulated, the aircraft has six state variables and three 
control inputs defined as 

x * [x y z V y 
u * [a $ T/m] T 


5 



6 


4 

H 






It 



K r 


5T 


Figure I.- Schematic Diagram of Aircraft Landing 
Simulation. 
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The state vector x is made up of the variables x, y, and z 
describing the aircraft's position; V, its velocity; y, its 
flight path angle; and j|>, its track angle. The control vector 
u is made up of the angle of attack, a; the roll angle, $; and the 
ratio of engine thrust to airc^ft mass, T/m. The two perturbation 
vectors 6x and 5u, shown in figure 1, are defined as 

5x = X - 

n 

<$u = u - u 
n 

However, in practice, 6u is computed as: 

H = KSx (2) 

where K is a 3x6 time varying gain matrix. This gai. matrix is the 
solution of the regulator problem which results from linearizing 
equation 1 about the nominal flight path and control vector, and it 
is precomputed for a particular nominal flight path. This 
linearization yields: 

fix - A(t)6x + B(t)6u * [A(t)+B(t)K(t)]6x (3) 

where A(t ) = t* evaluated at x = x and u = w and B(t) * r* 

3x n a 3u 

evaluates at x = x and u * u . A is a 6x6 matrix and B is a 6x3 
n n 

matrix. A statement of the regulator problem is then to find K in 
equations 2 and 3 which minimizes deviations of the states from 
their nominal values by using controls which stay within acceptable 
limits. More formally, the problem is to determine K which 
minimizes the quadratic performance index J given by 

J * — 6x T (t f )M6x(t f ) + j j f [6x(t) T Q6x(t) + 6u(t) T R5u(t)]dt (4) 

'o 

subject to the constraints of equation (3). M, Q, and R are weighting 


matrices. 
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Equations of Motion 

The inertial axes system (x, y, z) has its origin at the 
oesired touchdown point on the runway. The x axis is parallel to 
the runway and is positive in the direction of landing. The y axis 
is horizontal, perpindicular to the x axis and is positive to the 
right. The vertical z axis is perpindicular to both x and y, and 
is positive downward. The aircraft is treated as a point mass, 
subject to the forces of gravity, engine thrust, and aerodynamic 
lift and drag. The rotating axis system associated with the aircraft 
is chosen such that its x axis is alvays alined with the aircraft's 
velocity vector. Two Euler angle transformations relate the 
inertial axis system to the rotating axis system. They are 


’ COS 4* 

sin ip 

0 

-sin tj) 

cos 

0 

0 

0 

1 


cos y 

Ty * 0 

-sin y 


0 sin y1 

1 ° 

0 cos Y J 


The-se rotations are shown in Figure 2. This axes system does not 
roll with ti aircraft. Therefore, the y body axis remains 
horizc-'tta and the z body axis stays in a vertical plane with the 
velocity vector. The rotation vector resolved along the rotating 
axes system is 

w * [$ sin Y, -y , cos y] 
and the velocity vector is 

V « [V 0 of 
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The forces acting on the aircraft resolved along the rotating axes are 


Drag 

= [-D 0 

o] T 

Thrust 

= [Tcosa 

0 -Tsina] 

Gravity 

= T T. [0 
Y ♦ 

0 mg] T 

Lift 

= T, [0 0 

<J> 

-l] t 

where T^ - 

1 0 

0 

0 cos <}> 

-sin <J) 


0 sin <J> 

cos 4> 


The lift force is perpendicular to the velocity and rolls with the 
aircraft- The transformation resolves this force along the y and 
z rotating axes. The angle of attack, a* is defined as the angle be- 
tween the velocity vector and the x body axis of the aircraft. The 
dynamics of the aircraft are then described by equating the change in 
linear momentum to the sum of the applied forces 

m(V+ojxV) = Drag + Thrust + Gravity + Lift 


V TO " -D* " Tcosal rmg sin y 

m 0 + m Ivipcosy =0+0 + 0 

P, LVy OJ [-TsinaJ mg cos y 


* 0 * 
Lsin<J> 
-Lcos<$. 


To facilitate integration of these equations, they are written with 
the derivative of the states V, and y on the left hand side. 
These three equations sire combined with three kinematic relations 
which resolve the velocity along the x, y, and z inertial axis 


[T T. r 1 0 

y v 


The six equations of state are then 
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x * V cos y cos 
y * V cos y sin 4* 
z = V sin y 

(5) 

V = -D/m + (T/o)cosa + g sin Y 

Y = -(T/mV)sinot + (g/V)cosY-(L/mV)cos<J> 

\p = Lsin^/(mVcosY) 

Equation 5 is the expanded form of equation 1. 

Aerodynamics and Aircraft Parameters 

The aircraft simulated in this study is a Boeing 737-100 which 
is a small two-engined transport airplane designed primarily to 
operate from short runways over relatively short distances. It is 
described in references 9 and 10. The equations for the lift 
and drag forces are 


where C T and are the 

Jyl D 

S is the wing reference 
q are defined as 


L 2 q S C 

L 

D * q S C D 

coefficients of lift and drag, respectively; 
area; and q is dynamic pressure. C T , C_, and 

li .u 


C L = V a -°o>I30 

VS * " [C L (a - a o ) ilo l2 (6 > 

o a 

and q = 1/2 pV 2 

The assumption of a parabolic drag polar where the induced drag is a 
quadratic function of the lift coefficient is not restrictive for the range 
of speed and angle of attack employed in this study. 
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Air density, p, was approximated by the following equation taken from 
reference 8. 

p = 1.22[l.-(2.257xl0' 5 )h] 1 * -255 (7) 

This proved to be a very good approximation to the standard 

atmospheric density of reference 7* The airplane was configured with 

the landing gear down and the flaps deployed at 1*0 degrees. The 

parameters of equation 6 (C , a . C , and p) were chosen to fit 

Li 0 D 
a o 

data on and provided by Boeing for this configuration. The 

parameters of the aircraft are listed in Table 1. 

Nominal Flight Paths 

The nominal path can be chosen quite arbitrarily nen applying 
the linear regulator method to the aircraft landing problem. However, 
practical considerations lead one into choosing a particular class of 
flight paths. The most common approach used today on commercial 
transports is the 3° nondecelerating straight-in approach. That is, 
the pilot lines up with the runway and then descends at constant velocity 
along an equilibrium path which intersects the horizontal at an angle 
of 3°* The control inputs are nearly constant varying only to offset 

atmospheric disturbances. A mathematical description of this path is 

• 

T 

x = [V cosy 0. V sinY 0 0 0] 

o 'o o ’o 

x(0) * [x 0 2 V y n 0] T 

o o o o 

At approximately 12 meters altitude, the aircraft leaves this quasi- 
steady condition and flares in order to touch down more smoothly. This 
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SYMBOL 



n 


VALUE 

7.162 
- 10 . k° 
.157 
.031** 
91. OU M 2 
40,823. Kg. 


PARAMETER DEFINITIOH 
Lift curve slope 
Angle of attack for aero lift 
Drag coefficient for zero lift 
Efficiency factor 
Aerodynamic reference area 
Mass of vehicle 


TABLE 1 - Characteristics of the 737-100 (Gear Down-Flaps at U0°) 
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path calls for constant control inputs and because of its simplicity, 
it requires very little storage. If the nominal path had been chosen 
arbitrarily, the six states and three controls would have had to 
have been stored for many points along the path. 

In order to relieve airport noise problems and to decrease 
traffic congestion in the neighborhood of the airport, various other 
approaches have been proposed and are being tested. These include 
the steeper 6° nondecelerating straight-in approach and a combination 
of the 6° and 3° approach. The method prop'" jd in this paper is a 
quasi-steady helical approach. All other quasi-steady motions of an 
aircraft are special cases of this generalization. That is, the 
nominal path is descent along a helix from some initial altitude. 

The centerline of the runway is tangent to the helix and at some 
specified lower altitude, the aircraft follows this tangent into a 
straight-in approach. By making the assumption that atmospheric 
density is constant over the range of altitude considered, this path 
calls for fixed controls during the helical descent, and fixed but 
different control inputs during the straight-in portion. The 
variation in atmospheric density in this range of altitude is about 
four percent. The assumption of constant density is made during 
calculations of nominal path, nominal control inputs, and feedback 
gains. During the simulated flights, using the closed loop control 
system, atmospheric density is allowed to vary according to equation 
7. The small error introduced by this assumption was nulled out by 
the control system as though it were an external disturbance. A more 
formal description of the nominal path is 
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V cosy CO sill 

O 0 


“3c ~ 
0 

V o cosy o sin^ 


y c 

V siny 
o ’o 

0 

x(0) * 

z 

o 

V 

n 


0 

u 

ip 

y o 


y o 

\p 



L°J 


for the altitude range 


(8) 


-*SI i h i - z o 

where "Zgj is the altitude at which the aircraft rolls into the 

straight-in approach. The parameters z , z , V , y , and $ are 

O bl o o o 

chosen with some degree of latitude and they essentially determine 
the nominal flight path. The last three parameters, x q , y Q , and 
are functions of the first five and are chosen such that x, y, and ij> 
have the desired values when z equals z . The subscript SI will be 
used to specify the value of a state when z equals Zg^. Because the 
origin of the axes system is at the touchdown point on the runway, 
the desired values of ^gj. and yg^ are zero. In order to determine 
Xgj, it is necessary to look at the nominal path after it starts the 
straight-in approach. Its description is 

— Ip 

x * [V cosy 0 V siny 0 0 0] 

o o o o 

T <91 

x(t SI> = tx SI 0 Z SI V o t> 0] 


During this part of the flight x and z are constant and consequently, 


dx dx / dz 

T~ * 7T / 7T * eoty 

dz dt dt o 


is constant. In order to touch down at the origin. 
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X SI " Z SI dz 


( 10 ) 


X SI ■ 2 SI /tB "'<, 

Also, during the helical descent, the rate of turn, ip, and the rate 
of descent, z, are constant. In order to have the aircraft flying 
in the correct direction ( \p=0 ) when z equals Zg^, it is necessary 


that 


= 27T ^< 


dz 


2 sr z o 


and since 


it follows that 


dti> d\l> dz 
dt = dz dt 



( 11 ) 


( 12 ) 


The axis of the helix is vertical at 

(x,y) = ( x si» r ) 

where ;he ~adius r is given by 

r » l./(tanY • •jjj*) (13) 

o dz 

and through geometrical considerations 

X o * X SI + r 3in (lU) 

y Q * r(l-cos iji ) (15) 

A short example may clarify these computations. Let the landing 

direction be north, (\J>*0), and the aircraft's original heading be 

east, (<l» 0 *|’)« The velocity Is 62 meters/second and the flight path 
angle is 3°, (^q)* Initial altitude is U50 m. The aircraft descends 
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along a helical path to 50 m. altitude and then down a 

straight path to the runway. In this example, 

z = -1+50 m. , z _ = -50 m. , V = 62 m./sec 
o oi o 

Y = 22L_ W. , 2L 
y o 180 ’ ^o 2 

X SI = z SI /tan y o = - 50/tan( lfe ) * ~ 95k '*- 


difr _ o-n-ip _ 3tt 3tt = 

dz ” — “ 2(-50+U50) " 800 

z sr 7 „ 


.012 rad/m. 


^o = if It = { ioo )(62)sin Ho 2 * 038 rad/sec 

r = 1 /(tany o # ^-) * 800/(3ir*tan^Q) = l 620 .m. 


x = x OT + r sin to = -95 1 * + (1620;- in(-^—) = 666. m. 

O ox O c. 


y Q = r(l-cos(|» o ) = r - l620.m. 



Likewise, integration of equation (9^ yields the state trajectories 
along the straight portion of the flight path 
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Hie 


v 


V cosy (t-t ) 
o o SI 

0 


0 

Z SI 


V o s iny o (t-t S i) 

V 

o 

♦ 

0 

Y o 


0 

0 


0 

. 


L J 

r the 

visgs 

level maneuver. 


( 17 ) 


by 


t 


SI 



-z )/ u siny 
o o o 


(18) 


In sunenary, the need to integrate and store the equations of 
action for the nominal path has been eliminated. Given a particular 
input parameter set, [z q , z si , V q , y q , * 1, equations (16) and (17) 
yield the value of the states at ary time. Five secondary parameters, 
^ X SI* *o’ X o* y o* t si^* necessary to evaluate equations (l6) and 
(17) are given as functions of the input parameters by equations (10) 
through (15) and equation (18). 


nominal Control Inputs 

After specifying this particular class of nominal flight paths, 
it is necessary to make their equations of motion vhich are given by 
equations (8) and (9)» agree with the equations of motion for an 
arbitrary path which are given by equation (5). This is accomplished 
by choosing the correct open loop c .irol vector, u. Because the 



19 


path is quas i -steady , that is, 
V * y = 0 

$ = 0 



u is a constant vector on the spiral and a second constant on the 
straight-in portion. In order to make equation 15) agree with 
equation (8), the following equations must be satisfied. 


D H /m - (T/m)cosa - g sin y = 0 

n O 

(T/mV o )sina - (g/V^cosy^ + (L^/mY )cos+ = 0 (19) 

$> o - LjjSin$/(mV o cosY o ) = 0 


In these equations, 

hi = ^ S C L 

D H * 5 C D 

where 

- 1/2 »H V o 2 

and and are defined in equation (6). Using the assumption of 
constant air density discussed earlier, an intermediate value of 
altitude, h^, is chosen, and is computed using equation (?). 

That is, 

P N = 1.22[l.-(2.257xl0“ 5 )h H l 1 ** 255 


The open loop control, u, 

u = [a, ♦, T/m] 

which is the solution of equation (19)» is not unique and these 
equations have no closed form solution. Equation (19) is of 
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the form 


fi(u) * 0 

and it vas solved by the Newton-Raphson iterative method described 
in reference 11. This method usually requires the Jacobian matrix 
which is given by 


3D K 

^cT* + Tsina 3/° 



where 


3l n 


• (T/mV )cosoH-r-7Cos<J)/(inV ) 
o da o 


i 


L3T sin * /( mV o COSY o) 
3a 


0 

~L R sin<^/(mV o ) 
-L^cosfj)/ (mV Q cosY o ) 


— cosa 
sina/V 

o 


0 




3D„ » ? 

3a" = 2q N S ( llo ) C L a n(a “°o ) 

The possibility always exists that this algorithm will converge 
to an undesirable value of u. Physically, the smaller value of angle of 
attack is desired. This corresponds to the low side of the lift/drag 
curve. Therefore, the solution must be examined subjectively before 
being accepted. A practical way to avoid this problem is to start the 
iterative procedure close to the desired final solution. A suggested 
method of doing this is to make the assumption that engine thrust is 
alined with the velocity vector of the aircraft. W.'‘h this assumption, 
which is in error by the size of the angle, alpha, equation (19) becomes 

D N /m - T/m - g sinY 0 = 0 

-(g/ v 0 )c°sY o + (Ljj/mV^cos^ * 0 

^» o - L N sin<j>/(mV o cosY 0 ) * 0 
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These equations have the following closed fora solution which can 
he used as a starting point in the Hewton-Raphson algor it ha. The 
solution is 

♦ * Tan-^/g) 

a * o q ♦ ( 180/ 7T ) (BgcosY 0 ) / ( q^S Cj^cos^) ( 20 ) 

T/m * D_/m - gsiny 
« o 


In order to solve for the open loop controls for the straight-in 
portion of the nominal flight path, the following equations must he 
satisfied. They result from making equation (5) equivalent to 
equation (9)* 

D /m - (T/m)cosa - gsiny * 0 

N 0 (21) 

(T/mV o )sina - (g/V Q )cosY 0 + (V“V = 0 

Again, the Newton-Raphson algorithm is used to solve the two equations 
for a and T/m. The roll angle, +, of course, is zero. Using the same 
assumption that thrust and velocity are codirect ional , the following 


equations can he used to start the iterative procedure . 

a ■ lfiO.mgcosy /(irqS C. ) ♦ a 
o o 

T/m * Djj/m - gsiny o 


Notice that D R in equation 20 will have a different maerical value 
than Djj in equation 21. These equations used for starting the Hewton- 
Raphson iteration procedure yield values which are very dose to the 


final values. 
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Feedback Control Gains 

The purpose of this section is to present a develo pmen t of the 
solution of the linear regulator problem. References 1 through 5 
all concern themselves » in part* with this problem. A statement of 
the problem is given by equations (l), (3)» and (U), which are 
repeated hero for convenience. 


x - f(x,u,t) 


( 1 ) 


d:; = A(t)5x + B(t)6u = [A+BK]6x 


J = |dx T (t f )M6x(t f ) + 



6x^(t)Q6x(t 


+ 6u (t)R6u(t)dt 


(3) 

( 1 *) 


A(t) and B(t) are defined as 


A = 


3f 


B = 


3f 

3u 


M, Q, and R are weighting matrices which must be selected by the 
control system designer. The problem is to find the matrix K in 
equation (3) which minimizes the quadratic performance index of 
equation (U). Following the techniques of references 1 and 2, the 
Hamiltonian of the system is 

H = ^ 5x^Q5x ♦ j 6u*R6u ♦ (A6x)^p ♦ (B6u)**p 

where p is the vector of costates, p is the solution of the 
differential equation 

5 ’ ‘ 8Sr ' ^ <23) 


Since the control is optimal, it follows that 
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or 


lu *n_ 

sR5u>B ? s ° 


Solving this equation for "Su, ve get 


6u * -R _1 B T p (2U) 

Obviously, the control -eight ing matrix P. must be nonsingular. 
Combining equation (2kj with equation (3) yields 

6x - A6x - BR (25) 

Then, combining this equation with equation ( 23 ) » v e have the 
following system of equations. 


rfi 


-1 T~ 



! fix 


A -BR B 


fix 

i 

i 2. 


rp 



L P_ 

1 

1 

_-Q -A 

i 

- p. 


(26) 


Applying the transversality conditions, ve get the costate vector 
p(t) at the final time, t f . 

- H!t^ *®V <27) 

The system represented in equation (26) has tvelve states. Since 
there are six initial conditions on fix, and six final boundary values 
on p, this system has a unique solution. 

If equation ( 26 ) is rewritten as 


V 


6x 

• 

= Z 


_ p . 


_ P _ 


and if ve assume that Z is constant over a small interval of time, 
[t,t+h], then ve can write the transition equation of the system as 

^Reference 1 shows that M, Q, and R must be symetric; that M and Q must 
be at least positive semidefinite } and that R must be positive definite. 



2h 


where 



(28) 


(29) 


A further assumption is made that the costates are linearly related 
to the states; that is, 

p(t) = S 1 5x(t) (30) 


This assumption is proven to be true in reference 1. A comparison 
of equation (27) and (3c) shows that 

VV = M (31) 
Combining equations (28), ( 29 ), and ( 30 ) yields a transition 
equation for 


6x(t+h) 

S^(t+h )?x(t+h) 


fi(t,h) 


Sx(t) 

S 1 (t)?x(t) 


s 1 (t+h)[r. 11 +n 12 s 1 (t)]?x(t) * [ Wl ]©t) 


or 

S 1 (t+h) = [n 21 +fi 22 S 1 (t)][Q ;L 1 +n i 2 S 1 (t)]*' ;L (32) 

Reference 1 proves that this inverse does exist. 

Using the difference equation ( 32 ) vith the linal value of the 
matrix given by equation ( 31 ), we can proceed backward in time and 
obtain for all te[t o ,t f ). At each point in time, it is necessary to 
compute ft(t,h) as a function of A(t), B(t), R, and Q, using equations 
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(26) through (29). The gain matrix K(t) is obtained by combining 
equations (30) and (2U). 

6u * -R _1 B T S 1 6x(t ) 
or 

K(t) * -R' 1 B T S 1 (33) 


Implementation 

The landing approach guidance scheme proposed in this paper 
would be implemented in the following manner. 

1. For a particular nominal flight path, the control system 
designer, after extensive analysis and testing, chooses weighting 
matrices M, Q, and R for the performance index in equation (U). 

2. The state trajectories on the nominal flight path are then 
computed using equations (10) through (lQ). 

3. The nominal open loop control inputs are then computed using 
the Newton-Raphson algorithm to solve equations (19) and (21). This 
iterative procedure is started with the values given in equations (20) 
and (22), respectively. 

4. The Jacobian matrices A(t) and B(t) used in equation (3) are 
computed according to equations (A-l) and (A-2) in appendix A. Equation 
(3) represents the linearization of equation (l) about the nominal 
trajectory and control. 

5. The feedback gain matrix K is computed as a function of time 
using equations ( 26 ) through ( 33 ) * 

6. In the computer aboard the aircraft are stored the time his- 
tories of the nominal flight path, the nominal control inputs, and the 
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feedback gains. Alternatively, the algorithm, equations (l6) and (17), 
which computes the state trajectories is stored. 

7. The aircraft, or simulated aircraft, is flown into a window; 
i.e., a region of state space which is 'close' to the initial states 
of the nominal flight path. 

8. When the plane enters the vindow, it is switched on to the 
automatic landing system. It is assumed that the aircraft is 
receiving the necessary position data from the ground based airport 
landing system and receives control surface positions from transducers. 
This is considered time zero with respect to the nominal path. 

9- The control system smoothly nulls the state errors and brings 
the aircraft safely onto the nominal flight path well before reaching 


the decision altitude. 



RESULTS AND DISCUSSION 


Flight Paths 

In this study, the control system was tested along seven 
different nominal flight paths. The paths differed in initial 
altitude, initial heading angle, and in angle of descent. Flight 
paths 1, 2, and 3 are depicted in Figure 3- They all have a three 
degree angle of descent (flight path angle) and are initially headed 
north, east, ana south, respectively. Flight path 1 is the 
easiest path to fly because its initial altitude is the highest and 
there is more time to null out initial errors in the state variables. 
Conversely, flight path 3 is the most difficult since its initial 
altitude is 150 meters below flight path 2 ana 300 r.eters below 
flight path 1. The aircraft makes a descending turn to the right 
of 360, 270, and 180 degrees for flight paths 1, 2, and 3, respectively 
at which point it is headed north toward the runway. At this position 
(roll out) on the path, the aircraft is at 100 meters altitude and it 
rolls to a wing level attitude. It then descends in a straight line 
toward the runway. The flare maneuver would be executed at twelve 
meter altitude for the purpose of softening the impact at touchdown. 
Because this simulation did not include ground effects, it was 
terminated just above the flare altitude and off nominal errors 
recorded at that point. 

Flight paths 1A, 2A, and 3A had the same initial altitude and 
heading angle respectively as flight paths 1, 2, and 3. However, 
they had a descent angle of six degrees rather than three. These 
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flight paths are more difficult to fly than the first three. The rate 
of descent and roll angle are about doubled and the turning radius 
and flight time are about halved. The rollout altitude was maintained 
at 100 meters. These flight paths require more severe control inputs 
than the first three and they have less time to null off-nominal 
errors. 

Flight path 3B is a variation of 3A. The initial altitude and 
rollout altitude were raised 100 meters. This had the effect of 
moving the helix back from the runway and about doubling the time 
spent on the straight-in portion of the path. Parameters of the 
seven flight paths are given in table 2. 

Design Considerations 

The control task presented here differs from conventional landing 
approach schemes in that errors are taken with respect to a moving 
point in state space rather than a fixed line. Thus, even if the 
aircraft were flying on the correct helical path • ith the correct 
attitude, the fact that it was late or early in time would indicate 
to the control system that there were state errors; in this context, 
the proposed guidance scheme can be considered U-D. At crowded 
airports, more precise contro 1 of time sequencing should prove 
advantageous. 

Because the feedback gain matrix K is computed automatically in 
a near optimum fashion, the designer's Job is to choose the weighting 
matrices M, Q, and R of equation 1* and then to evaluate the subsequent 
performance of the control system. In fact, this becomes an iterative 
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procedure where the designer is trying to find an optimum set of 
weighting parameters- Tnis task is not trivial, but by using 
diagonal weighting matrices, a functional relationship between their 
elements and the subsequent performance can be established. 

The control system was required to null out initial errors in 
the states, and to offset the effects of steady crosswinds. It was also 
desired that it not be sensitive to normal changes v the atmosphere 
and uncertainties in tne aircraft's aerodynamics. Wind gusts and 
sensor noise were not considered in this study. Large initial state 
errors tend to demand large rapid control motions. In order to avoid 
tius, state errors ar^ weighed less by reducing the elements of Q and 

control excursions are weighed more by increasing the elements of R. 

However, the penalty function does not include control rates and the 
R matrix influences tne time integral rather than the magnitude ^f 
control excursion squared. Therefore, it is not unusual tc j *oge 
control excursions for relatively short time with consequent high 

control rates. If the Q and R matrices are well chosen, the control 

system will null out the initial errors slowly in a manner which will 
cause no passenger discomfort. As the aircraft approaches the 
touchdown point, the influence of the M matrix tightens up the control 
and acceptable values of final state errors can be achieved. 

This approach worked fine until the aircraft was subjected to a 
crosswind. It then became apparent that control in crossvinds was a 
much more stringent requirement and thet a tighter control was 
necessary. In order to reconcile the two requirements: vide 


bandwidth control for crosswinds and narrow bandwidth control for 
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large initial errors, two approaches were taken. First, high gains 
were used, hut rate and position limits were placed on the controls. 
This lead to undesirable oscillatory behavior, and in some cases, 
instability. The rate limit on angle of attack seemed to be 
pa? icularly troublesome. The second approach was to make the 
elements of Q and R time varying. In this method, the feedback gains 
are low initially, while the state errors are large. While the large 
errors are being nulled, the feedback gains are gradually increased 
and the system is reasonably tight when it approaches touchdown . 

This second method works well except that roll rate became excessive 
when the aircraft rolled out of its turn. This maneuver occurs at 
100 meters altitude when high feedback gains are desirable. Thi3 
proolem was circumvented by imposing a rate limit on roll . The 
values of Q, R, and M which were finally settled upon are: 




33 



where the value of the elements are given fcy 


5x q = 

i^O 85 - t/t^. 

s = 

lU0.-131-t/t f 

6z q = 

3 i.- 26 t/t f 

5v q • 

70.-l*9.t/t f 

^ = 

3 - (tt/180 . ) 

% 

( 7 .- 6 .t/t f )( 7 r/i 80 ) 

6a B 

3+12t/t f 

6^ = 

(1.7ir/l80)(i+t/t f ) 

5(T/M) r = 

.17+.17t/t f 

6x m 

55- 

6y M 

9- 

II 

SS 

*o 

5- 

6v m = 

21. 

6y m 

2. (tt/180) 


tt /180 


These weights are by no means optimally chosen. They are arrived at 
experimentally by using the iterative process described above. 



There axe several constraints imposed upon the design. To 
avoid passenger discomfort, angular rates should not exceed ten 
degrees per second and by federal regulation, roll angle cannot 
exceed thirty degrees. Engine thrust can be throttled between 2800 
and 28,000 pounds which correspond to thrust /mass ratios of .3 and 
3- meters per second squared. Final errors in the states will be a 
function of how well the weighting matrices are chosen, the magnitude 
of the initial state errors, wind velocity and direction, and of the 
particular flight path chosen. After the wings level maneuver the 
aircraft is headed north at about 62 meters/second, and since the 
runways are normally much longer than necessary for this aircraft, 
an error in x of + 62 meters is not unreasonable. However, errors 
in y of greater than +. 10 meters are unacceptable. The fact that the 
runway is normally only 50 meters vide and the wing span is approxi- 
mately 30 meters make even that error uncomfortable even though some 
correction will be made during flare. Errors in velocity V and 
flight path angle y are best viewed in terms of sink rate Z 
(Z = V sin y). This error should be kept within .7 meters/second 
and again a negative error (slower descent) is preferred. Error in 
heading should be kept within five degrees. An error of that 
size can be nulled by a decrab maneuver. The most trouble- 
some errors in this study turned out to be in y during the six degree 
approaches . 

Representative Landing Approach 

Time histories from a representative flight are shown in figure U. 
The nominal curves axe shown as dotted lines and the simulated flight 
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&s a continuous curve. This was a six degree approach with nominal 
initial heading due east and initial altitude of 550 meters. The 
six states and three controls are shown, along with an overview of 
the x-y plane. The aircraft started wings level with a heading of 
80° (10° error). It was 100 meters higher than nominal and had 
over a 300 meter error in both x and y. The velocity was low by 
about twenty percent. The aircraft was near stall, and consequently, 
the angle of attack and throttle setting were set high for trim. 

It was subjected to a fifteen knot constant east wind. Immediately, 
the aircraft rolled with some overshoot to the 23° nominal. The 
angle of attack was decreased and the throttle was cut back. This 
made the aircraft simultaneously pitch over and increase speed 
while initiating its curved descending turn to the right. The roll- 
out at 100 meters altitude is apparent in the curves of y, psi, and 
phi. The entire approach lasted about 82 seconds and rollout 
occurred at around 70 seconds. The simulation was terminated when 
the plane reached an altitude of around 15 meters just above where 
the flare would be executed. Most of the notion is smooth but roll 
angle rate definately reaches its ten degree per second limit and the 
thrust to mass ratio stays on its lower limit of -3 (T - 2800 lbs.) 
for much of the run. The overshoot in roll angle (-10°) at rollout 
is undesirable, but is much less apparent in the three degree 
approaches. It could, most likely, be considerably reduced by making 
the nominal roll angle a continuous function. However, this would 
complicate (perhaps unjustifiably) the entire control system. At 
decision height, the aircraft was within one second of the nomin&x. 
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Final errors for this flight are -52m., 3m., -.6n-, 3 m./s. f -1°, 

.5°, and -.5m./s. for x, y, Z, V, y, \p, and Z, respectively. 

Representative Feedback Gains 

The elements of the feedback gain matrix, K, for flight path 
2A are shown in Figure 5* The magnitude of the gains, in general, 
increase with time due to the time variation of the weighting 
matrices Q and R. Some of the gains have a pronounced increase in 
the last several seconds of flight due tc the influence of the M 
matrix* This can be seen in several of the elements including K(l»l) y 
K(2,2), and K(3»3). During the straight-in portion of the approach, 
ten of the eighteen gains are zero or can be considered zero. Of 
course, the more of these elements that are zero or near zero, the 
easier it is to implement the control system. In order to determine 
the influence of a particular gain, one assumes a large reasonable 
error in the state and then calculates the resulting change in the 
control. If the change is insignificant , then the element can be 
considered zero. For example, an error of three meters in Z near 
the end of the flight would only change the thrust /mass ratio by . 06 . 
Whereas an error in x of 60 meters would cause the ratio to change by 
1.2. Since these are both large reasonable errors at the termination 
of flight, the element K(3,3) should be considered zero. Using this 
approach, it can be seen that after rollout, thrust/m&ss ratio is a 
function only of x and V, roll angle is a function of only y and \p 9 
and angle of attack is a function of x, V, Z, and y. This partial 
decoupling is useful to the designer in the iterative process of 
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choosing the weighting matrices. Applying the same reasoning, five 
elements of K can be consider el zero throughout the entire flight. 

They are K(2,3), K(2,M, K(2,5), K(3,5), and K(3,6). Each of the 
remaining 13 gains would have to stored in the on-board computer. 
Eleven points with linear interpolation could easily approximate any 
of these functions. Some of them could be fitted by low order 
polynomials in time which would require even less parameter storage. 

It appears that the entire gain matrix and nominal path .culd be 
stored in less than two hundi d words of storage. This number could 
probably be reiuced but it is important to realize that any para- 
meter change in the nominal state trajectory wilx affect these gains. 
Each airport using this system wotild probably have several different 
nominal approaches ai.a these would also probably differ among airports. 
To be practical then, the computer must have access to the correct set 
of parameters out of many possible alternative sets. Several 
possible ways of accomplishing this are: 

1) All necessary parameter sets are stored .n an on-board 
mass storage device. 

2) The particular set of parameters needed is calculated 
on the airborne compute/ prior to landing approach. 

3) The parameter set is transmitted from the airport to 
the on-board computer when the landing approach is 
assigned. 

k) The feedback gains could possibly be expressed as 
analytic functions cf time and vhe five primary 
parameters of the nominal flight path. 



The feedback gains of Figure 5 were calculated using the 
weighting matrices M, Q, and R as given above. The gains were used 
in the landing approach shown in Figure U. This choice of weighting 
matrices and resulting gains appears adequate for this flight path 
and control task. 

System Performance 

The control system was tested on each of the seven landing 
approaches described in Table 2. Results from the first six flight 
paths are given in Table 3 through Table 8. Flight path 3B results 
are given in Table 9* Each of these seven tables represent 25 
simulated landing approaches. On each approach, the aircraft was 
subjected to a constant wind of 15 knots magnitude from a direction 
^ . The direction changed for each approach and it was generated as 
a random variable with uniform distribution between -180 degrees and 
+180 degrees. Each flight path had its own random sequence of wind 
a rections with no particular direction being favored. Each approach 
was started with off nominal errors in all state variables. These 
initial errors were generated as random variables. They were nor- 
mally (Gaussian) distributed with zero mean. The standard deviations 
selected for the six state errors in x, y, Z, V, y* respectively, 
were 100m., 100m., 30m., 3m. /s., * \ and 3°. Each flight p tb had a 
distinct random sequence of sets of initial state errors. 

The reason for employing randomness in the control task was that 
the system showed that it could be tuned to a particular task. 
Originally, only 90° crosswinds were used because they were thought 


to be a worse case. Then it was discovered that winds from the aft 



quarter occasionally gave more trouble. Because this aircraft is 
relatively clean; i.e., high lift to drag ratio, these winds were par- 
ticularly difficult on the steep approaches. Also, certain combinations 
of initial errors were easier to correct than others • Therefore, a 
strong wind randomly directed was used and the aircraft initialized 
with large random off -nominal errors. The controls were initially set 
such that the wings were level ($ * 0) and the aircraft was trimmed. 

Table 3 presents the data for the 25 landing approaches along 
flight path one. The first six lines give the measured statistics 
of the initial state errors. They are the mean, standard deviation, 
and range of initial errors given for each otate. These statistics 
differ fror; their ideal values because they represent a subset of a 
very long pseudo random sequence which has the selected statistics. 

The next three lines of the table represent the initial off-nominal 
val.es of the controls. These are not really errors, but show a 
variation in initial trim values. Roll angle is a constant -12.2° 
since the approach always starts wings level rather than rolled over 
at 12.2°. The next section of the table shows how well the control 
system did its work. It gives the final off-nominal errors taken 
Just above the flare altitude. Statistics are given on the six 
state variables and sink rate, Z. As can be seen, the control system 
worked very well on flight path one. Poor performance usually makes 
itself apparent in final errors in y end Z. For these twenty-five 
runs, the approaching aircraft was never more than .7 meters from 
the center line of the runway and the final error in altitude was 
never greater than 1.1 meters. Since the nominal sink rate is 3.2U 
meters/second, the maximum error in final altitude corresponds to 
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Variable 

Mean 

Standard Deviation 

Range 



INITIAL OFF-NOMINAL ERRORS 

X (m) 

-6.9 

64.7 

(-146., 119.) 


-13.0 

84.5 

(-153- , 193.) 

Z (m) 

-4.1 

32.1 

(-85., 76.) 

v (m/s) 

.6 

2.8 

(-4.2, 5.2) 

Y (deg) 

.2 

.8 

(-1.3, 1-9) 

4 > (deg) 

■ 5 

3.1 

(-6.1, 8.6) 

a (deg) 

.1 

1.4 

(-2.1, 2.7) 

<J> (deg) 

-12.2 

0 


T/m (m/s ) 

-.05 

.Ik 

(-.36, .20) 


FINAL OFF-NOMINAL ERRORS (AT DECISION ALTITUDE) 

1 

x (m) 

-9-3 

7-8 

(-21.2, -1.0) 

y (m) 

.1* 

-2 

(.3, .7) 

Z (m) 

.0 

.6 

(-.6, 1.1) 

V (m/s) 

.2 

-5 

(-.3, -9) 

Y (deg) 

-.0 

.3 

(-.5, .M 

(deg) 

-.1 

.1 

(-.2, -.1) 

Z (m/s) 

-.01 

• 34 

(-.47, -44) 


EXTREME OFF-NOMINAL EXCURSIONS DURING EACH FLIGHT 

1 

a (deg) 

4-5 

• 5 

(3.3, 5-7) 

♦ (deg) 

-.1* 

2.3 

(-2.3, 3.9) 

T/m (m/s ) 

.08 

.18 

(-.35, .30) 

! 

EXTREME CONTROL RATE DURING EACH FLIGHT 

a (deg/s) 

.1* 

1.5 

(-2.7, 3.4) 

* (deg/s) 

7-9 

3.1 

(-6.7, 9-4) 

T/m (m/s^) 

! 

.00 

.11 

(-.21, .22) 


NUMBER OF UNACCEPTABLE FLIGHTS = 0 


TABLE 3.- Statistical Data From Flight Path 1. 

25 Landing Approaches, y* 3°» Z*3.24 m/s, V *15 knots, 
<J> w =U[-l80°, l80°] , t f =211.5s, h f =13.7m, x f *-262m. 






— 

Variable 

— 

Mean 

Standard Deviation 

Range 



INITIAL OFF-NOMINAL ERRORS 

x (m) 

8.4 

118.1 

(-292., 236.) 

y (m) 

10.5 

93-9 

(-199., 179- ) 

Z (m) 

2.3 

21.5 

(-27-, 48.) 

V (m/s) 

.2 

2.7 

(-4.8, 6.9) 

Y (deg) 

.2 

• 9 

(-1.4, 2 . 9 ) 

<i> (deg) 

-•9 

2.2 

(-5-7* 3-0) 

a (deg) 

.1 

1-3 

(-2.9, 2.8) 

<t> (deg) ? 

-12.2 

0 


T/m (m/s ) 

-.06 

.14 

(-1*6, .19) 


FINAL OFF-NOMINAL ERRORS (AT DECISION ALTITUDE) 

x (m) 

-8.9 

6.8 

(-22.9, --!*) 

y (m) 

• 5 

.2 

(.2, .8) 

Z (m) 

-.1 

.6 

(-.7, 1.2) 

V (m/s) 

.2 

• 3 

(-.3, .9) 

Y (deg) 

-.0 

.2 

(-5, A) 

^ (deg) 


.1 

(-.2, .0) 

Z (m/s) 

-.02 

-29 

-3* 

* 

GO 

• 

1 


EXTREME OFF-NOMINAL EXCURSIONS DURING EACH FLIGHT 

a (deg) 

2.9 

3-2 

(-3-7, 6.6) 

♦ (deg) 

4.2 

.1 

(U.l, U.U) 

T/m (m/s^) 

.06 

.17 

(-35, -29) 


EXTREME CONTROL RATE DURING EACH FLIGHT | 

! 

a (deg/s) 

•5 

2.0 

(-3.0, 4.8) 

6 (deg/s) 

8.8 

• 7 

(6.7, 9-4) 

T/m (m/s^) 

.01 

.15 

(-.29, .23) 


NUMBER OF UNACCEPTABLE FLIGHTS = 0 


TABLE 4.- Statistical Data From Flight Path 2. 

25 Landing Approaches, y= 3°» 4=3- 24 m/s, V w =15 knots, 
^“Ut-180 0 , l80°], t f *l65.s, h^=l4.6m, x f *-279*in 
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Variable 


x (m) 

y (m) 

Z (m) 

V (m/s) 

Y (deg) 
tj> (deg) 

a (deg) 

$ (deg) ? 
T/m (m/s ) 


x (m) 

y (m) 

Z (m) 

V (m/s) 

Y (deg) 
^ (deg) 
Z (m/s) 


a (deg) 

♦ (deg) 
T/m (m/s^) 


(deg/s) 

(deg/s) 

/m (m/s^) 


Mean 

Standard Deviation 

Range 


INITIAL OF NOMINAL ERRORS 

11.1 

95.2 

(-176., 171*. ) 

-12.8 

104.1 

(-209., 217.) 

-2.5 

30.8 

(-48., 73.) 

-.0 

2.1 

(-3.7, 5.2) 

-.1 

1.1 

(-1.8, 1.8) 

.1* 

2.5 

(-4.1, 6.2; 

-.1 

1.0 

(-2.4, 1.8) 

-12.2 

0 


-.02 

.20 

(-.35, .29) 

FINAL OFF-NOMINAL EF.RORS (AT DECISION ALTITUDE) 

-8.8 

9-3 

(- 26 . 0 , 2.2) 1 

.6 

.3 

(. 0 , 1 . 3 ) ' 

.1 

• 5 

(-5, 1 . 0 ) 

.1 

.6 

(-.7, 1.2) 

.0 

.2 

(-.2, .4) 

-.1 

.1 

(-.3, -. 0 ) 

.06 

.26 

(-25, .53) 

EXTREME OFF-NOMINAL EXCURSIONS DURING EACH FLIGHT 

2.8 

3.1 

(-3.8, 7-2) 

7-2 

.2 

(6-9, 7-5) 

.08 

.19 

(-20, .39) 

EXTREME CONTROL RATE DURING EACH 

FLIGHT 

-.2 

1.8 

(-3-5, 2.6) 

8.3 

• 9 

(6.4, 9-4) 

-.00 

.13 

(-.24, .20) 


NUMBER OF UNACCEPTABLE FLIGHTS = 0 


TABLE 5*- Statistical Data From Flight Path 3 • 

25 Landing Approaches, Y = 3°, Z=3-24 m/s, V =15 knots 
^ w =U[-l80°, 180° ], t f =ll8.5s, h f -15-5m, x f =-296m. 
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Variable 


Mean 


Standard Deviation 


Range 


INITIAL OFF-NOMINAL ERRORS 


X 

la) 

27 . u 

112.7 

(-219., 276. 

y 

(m) 

13-3 

106.0 

( - 22 *» . , 237- 

z 

(m) 

-7-7 

26.7 

(- 60 ., 1 * 0 ) 

V 

(m/s) 

-.1 

3.6 

(-7.2, 6.k) 

Y 

(deg) 

.1 

1.1 

(-1.5. 2.U) 

V 

(deg) 

-.3 

3-3 

(-6.8, 6.3) 

a 

(deg) 

-•5 

1.9 

(-3.1*. 3-7) 

♦ 

(deg) 

-23.3 

0 

(-.17, .17) 

T/m (m/s ) 

-.07 

.11 



FINAL 

0 

k 

1 
tr* 

1 

(AT DECISION ALTITUDE) 

X 

(m) 

-10.0 

9-U 

(-23.6, 8.3) 

y 

(m) 

8.1 

.7 

(6.6, 9.1) 

z 

(m) 

-6 

.6 

(-1.2, .7) 

V 

(m/s) 

.0 

.5 

(-1.1, .8) 

Y 

(deg) 

.1 

.2 

(-.2, .3) 

♦ 

(deg) 

-1.8 

.2 

(-2.0, -l.M 

Z 

(m/s) 

.11 

.25 

(-36, .U6) 


EXTREME OFF-NOMINAL EXCURSIONS DURING EACH FLIGHT 


c (deg) 

♦ (deg) 0 
T/m (m/s^) 


3.1 
19.3 
• 17 


3.8 (-3-7, 8.U) 

.1 (19-0, 19-5) 

•lfc (-.17, .35) 


EXTREME CONTROL RATE DURING EACH FLIGHT 


L 


(deg/s) 
(deg/s) 
/m (m/s ^) 


.1 

• 10.0 

.01 


2.1 (-3.U, 3.3) 

0 

.10 (-.18, .21) 


NUMBER OF UNACCEPTABLE FLIGHTS * 0 


TABLE 6.- Statistical Data From Flight Path 1A. 

25 Landing Approaches. y=6° , 1*6 . 28 m/s, V *15 knots 
4» w *U[-l80°, 180° ], t f =105.5s, h f =19 • 5m , x f *^l86m. 
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Variable 

Mean 

Standard Deviation 

Range 



INITIAL OFF-NOMINAL ERRORS 

x (m) 

- 9-7 

89.2 

(- 235 ., H * 7 .) 

i y v * a ) 

20.8 

75 - 1 * 

(- 11 * 3 -, 192 .) 

z ( a ) 

1.2 

26.2 

(- 1 * 1 ., 56 .) 

i V (m/s) 

• 3 

3.1* 

(- 6 . 9 , 6 . 2 ) 

| Y (deg) 

-.0 

.8 

(- 2 . 1 , 1 . 1 *) 

i ^ ( deg) 

-.3 

2.5 

(- 3 . 5 , 6 . 9 ) 

a (deg) 

-.9 

1.7 

(- 3 . 6 , 3 . 2 ) 

< J > (deg) _ 

- 23.3 

0 


T/m (m/s ) 

-.07 

.11 

(-. 17 , . 2 L ) 


FINAL OFF-NOMINAL ERRORS (AT DECISION ALTITUDE) 

x (m) 

- 1*.9 

9.0 

(- 25 - 1 *, 7 - 9 ) 

y (m) 

8 .L 

1.3 

( 6 . 5 > 10 . 6 ) 

Z (m) 

-. 1 * 

1.0 

(- 1 . 7 , 1 . 3 ) 

V (m/s) 

-.3 

.5 

(- 1 . 1 , - 9 ) 

Y (deg) 

.0 

• 3 

(- 1 *, . 1 *) 

^ (deg) 

- 1.7 

.3 

(- 2 . 1 , - 1 . 3 ) 

Z (m/s) 

.01 

• 25 

(-. 38 , . 33 ) 


EXTREME OFF-NOMINAL EXCURSIONS DURING EACH FLIGHT 

a (deg) 

3.0 

3-5 

( - 1 * . 3 , 6 . 5 ) 

4> (deg) 

13-3 

.3 

( 12 . 8 , ll ». l ) 

T/m (m/s^) 

.09 

• 15 

(-. 16 , . 26 ) 


EXTREME CONTROL RATE DURING EACH FLIGHT 

a (deg/s) 

-.1 

2.0 

(- 3 - 3 , 3 - 9 ) 

A (deg/s) 

- 10.0 

0 . 


T/m (m/s^) 

-.01 

.11 

(-. 17 , . 21 ) 


NUMBER OF UNACCEPTABLE FLIGHTS = 5 


TABLE 7-~ Statistical Data From Flight Path 2A. 

25 Landing Approaches, y =6°, Z=6.28 m/s, V *15 knots 
<J» w =U[-l80°, 180° ] , t f *82.5s, h f *15.3m, x f *-Yl*6.m. 
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Variable 

Mean 

Standard Deviation 

Range 



INITIAL OFF-NOMINAL ERRORS 

(a) 

-35-6 

88.L 

(-225., 129.) 

(») 

-22.2 

89.5 

(- 266 ., 81.) 

(m) 

-5.2 

36.3 

(- 61 ., 5 1 *-) 


.1 

2.5 

(-6.0, 1* .2) 

(deg) 

.1 

.9 

(-1.3, 2.3) 

(deg) 

-•3 

2.8 

(-5- 1 *, l*-7 ) 

(deg) 

-1.1 

1.3 

(-3.0, 2.U ) 

(deg) 

-23.3 

0* 


7m (m/s ) 

-.07 

.09 

(-.17, -09; 


FINAL OFF-NOMINAL ERRORS (AT DECISION ALTITUDE) 

(m) 

-21.5 

11*. 7 

(-51.7, 2.1*) 

(m) 

11.? 

5.5 

(-6.1, 20.3) 

(m) 

-.5 

• 7 

(-1.8, 1.2) 

(m/s) 

1.1 

1.0 

(-.8, 2.8) 

(deg) 

.1 

.2 

(-.1*, .3) 

(deg) 

-2.1 

1.1 

(-I 4 . 0 , l.l*) 

(m/s) 

• 17 

• 25 

(-.22, .57) 


EXTREME OFF-NOMINAL EXCURSIONS DURING EACH FLIGHT 

» (deg) 

• 5 

1* .2 

-* 

** 

CVJ 

• 

-=* 

1 

► (deg) 

21.7 

1.0 

( 19 . 7 , 23 . 1 *) 

P/m (m/s^) 

• 19 

.21 

(-.17, .55) 


EXTREME CONTROL RATE DURING EACH FLIGHT 

f (deg/s) 

1.1 

1.8 

(-2.0, 5.1) 

f (deg/s) 

-10.0 

0. 


L'/m (m/s^) 

.02 

.09 

(-.15, .15) 


NUMBER OF UNACCEPTABLE FLIGHTS = lU 


TABLE 8.- Statistical Data From Flight Path 3A. 

25 Landing Approaches. y=6°, Z=6.28 m/s, V *15 knots 
4» w =U[-l80°, 180], t =59. Os, h «15.n, x f =-lU0m. 









about one third of a second in time. 'In is same observation can be 
made with regard to x since nominal x is approximately 62 meters/ 
second and the maximum final error in x for the twenty-five runs 
was 21.2 meters. The third section of the table indicates bow much 
control was needed to fly the aircraft. For each of the twenty-five 
approaches, the largest off-nominal control excursion is identified 
whether it be positive or negative. The initial transient of each 
run was not included since the initial off-nominal values of at 
least roll angle would always be extreme for the flight and it 
would mask the desired information. So, for example, the average 
extreme for angle of attack was 1».5 degrees greater than the nominal 
and the extreme varied over the twenty-five flights from 3-3 degrees 
to 5-7 degrees above nominal. The next section of the table gives 
information on extreme control rates for each run. These are 
important mainly in how they affect passenger comfort. The 
statistics given are on the extremes of each flight including the 
initial transient. Finally, the table gives the number of flights 
which v’re unacceptable because one or more of the constraints, 
discr.S'Jed in the section on design consideration, was violated. 

As car. be seen from examining the tables, the control system 
was luite effective on flight paths 1, 2, and 3> The largest final 
errors for these approaches was 1.3 meters in y and 1.2 meters 
in Z. There were no unacceptable flights and both control 
excursion and control rates were moderate. 

The next three tables, 6 through 8, show the results for the 
six degree approaches. As can be seen, the control system did not 
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do as well on these approaches. Final errors in y and i|/ were larger, 
and roll angle tended to have larger off-nominal excursions. 

Although there were no unacceptable flights on path 1A, the mean 
final error in y was a disturbing 8.1 meters. Because the aircraft 
did not execute the wings level maneuver fast enough, it tended to 
go off to the right (east), and then it would roll to the left in 
order to get back on the centerline. On these six degree approaches, 
there was only about 13 seconds from initiation of rollout to the 
conclusion of the simulation. Since the variation of the error is 
small (.7), performance on this flight path could probably be 
improved by starting the rollout a few seconds early. This is, in 
fact, what pilots have been observed to do when flying this approach. 
On flight path 2A, performance is slightly worse and five of these 
had final errors in y which were greater than ten meters and 
consequently were unacceptable. The variation in final y error was 
again small, indicating that a slight adjustment in the flight path 
could possibly bring the errors within plus or minus 2.5 meters. 
Finally, on flight path 3A, the control system proved to be 
completely inadequate with lU of the 25 flights having final errors 
in y which were greater than ten meters. It is possible that another 
choice of weights, specifically 6y n and 6y„, might improve performance 
on the six degree flight paths. However f the principle problem seems 
to be one of time. Flight time decreases monotonically from 212 
seconds on flight path 1 to around 60 seconds on flight path 3A. It 
appears that the six degree approaches do not have enough time to 
null out initial errors and the adverse effect of crosswinds. In 
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order to support this premise, flight path 3A was modified such that 
initial altitude was 500 meters instead of 400, and rollout was 
accomplished at 200 meters instead of 100. This new flight path is 
labeled 3B. Its parameters are given in Table 2 and its performance 
is reported in Table 9. In order to facilitate comparison of Table 
9 and Table 8, the exact same sequences of initial errors and wind 
directions were used. The flight time for 3B was about 15 seconds 
longer than 3A, and all of + his was added to the straight-in portion 
of the flight . The marked improvement in performance of 3B over 3A 
is apparent when Table 9 is compared with Table 8. On 3B, no 
approach exceeded the allowable limit of ten meters final error in y; 
whereas on 3A, ll* approaches exceeded it. Since this improvement 
seems to be caused by the increased flying time, it may follow that 
doubling the initial nominal altitude end rollout altitude on the 
six degree approaches would lead to performance comparable to that 
experienced on the three degree approaches. 

Finally, the control system was tested for sensitivity to 
uncertainties in the aerodynamics of the aircraft and for normal 
variation in the atmosphere. A particular landing approach was used. 
Initial errors and wind direction were not varied. Twenty-five 
simulation runs were made, each with a different atmospheric density 
function and each with a different set of aerodynamic parameters . The 
results of these flights are given in Table 10. All six aircraft 
parameters were allowed to vary randomly. The values given in Table 1 
were used as the mean and four percent of those values were used as 


the standard deviation. The distribution was normal. This resulted 
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Variable 

Mean 

Standard Deviation 

Range 



INITIAL OFF-NOMINAL ERRORS 

x (m) 

-35.6 

88. U 

(-225., 129.) 

y U) 

-22.2 

89-5 

( - 266 . , 8l. ) 

Z (m) 

-5-2 

36.3 

(- 61 ., 54.) 

V (m/s) 

.1 

2.5 

(-6.0, 4.2) 

y ' 'eg) 

.1 

.9 

(-1.3, 2.3) 

vdeg) 

-.3 

2.8 

(-5-1*, 4.7) 

a (deg) 

-.9 

1.3 

(-2.9, 2.6) 

4> (deg) 

-23-3 

0. 


T/m (m/s Z ) 

- "7 

.09 

(-.17, .09) 


FINAL OFF-N ' L ERRORS (AT DECISION ALTITUDE) 

x (m) 

- 16.2 

ll». 2 

(-42.4, 13.2) 

y (m) 

4.8 

2.1 

(-.6, 8.5) 

Z (m) 

-7 

• 5 

(-1.1*, .6) 

V (m/s) 

.8 

.6 

(-.!», 2.1) 

y (deg) 

.0 

.2 

(-.6, .3) 

^ (deg) 

-.8 

.3 

(-1.1*, .1) 

Z (m/s) 

.12 

.3 

(-.54, .50) 


EXTREME OFF-NOMIKAL EXCURSIONS DURING EACH FLIGHT 

a (deg) 

.9 

li.l 

(-4.1, 7.2) 

♦ (deg) 

21.6 

.8 

(19.6, 23-1) 

T/m (m/s*) 

.18 

• 19 

(-.17, .1*6) 

! 

EXTREME CONTROL RATE DURING EACH FLIGHT 

4 (deg/s) 

1.1 

1.8 

(-2.0, 5.0) 

♦ (deg/s) 
T/m (m/sd) 

-10.0 

.0 

0. 

.1 

(-.15, .15) 


HUMBER OF UNACCEPTABLE PLIGHTS * 0 


TABLE 9-- Statistical Data From Flight Path 3B. 

25 Landing Approaches. y*6 ° , 4*6.28 m/s, V *15 knots 
*U[-l80°, l80°] , t *75.0s, h f «ll».7n, x *-l35*. 
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Variable 

Mean 

Standard Deviation 

Range 



1 

INITIAL OFF-NOMINAL ERRORS j 

x (m) 

186.0 



y (m) 

157-9 



Z (m) 

-2k. 0 



V (m/s) 

l.k 



Y Ueg) 

-.2 



ij> (deg) 

2.1 



a (deg) 

-.6 



4> (deg) ? 

-12.2 



T/m (m/ s ) 

.01 




FINAL OFF-NOMINAL ERRORS (AT DECISION ALTITUDE) 

X (m) 

-2.0 

8.4 

(-23.7, 12.4) 

y (m) 

.3 

.1 

(.2, .5) 1 

Z (m) 

1.3 

.8 

(-.2, 3.0) 

V (m/s) 

-.2 

• 5 

1 

M 

• 

M 

• 

O 

Y (deg) 

.k 

.1 

(-.0, .5) 

(deg) 

-.1 

.0 

(-.1, -.1) 

Z (m/s) 

i 

•43 

• 1 

(.05, .51) 


EXTREME OFF-NOMINAL EXCURSIONS DURING EACH FLIGHT 

) 

4.7 

.8 

(3-7, 7-1) 

) 0 

4.3 

.0 

(4.2, 4.3) 

i/s 2 ) 

-.lk 

• 15 

(-.32, .34) 


EXTREME CONTROL RATE DURING EACH FLIGHT 

4 (deg/s) 

-1.0 

• 3 

(-1.3, .5) 

* (deg/s) 
T/m (m/s^) 

8.9 

-.16 

.0 

.01 

(8.9, 8.9) 
(-.17, 0.15) 


NUMBER OF UNACCEPTABLE FLIGHTS * 0 


TABLE 10.- Statistical Data From Flight Path 2 With Variations in the 
Aerodynamic Parameters. 

25 Landing Approaches. Y*3°, Z*3.24 m/s, V w »15 knots, 
^■170°, t «l65., h *l4.6m, Xf*-279«». Random variation 
of Approximately ♦ 10# .n Aircraft Parameters and ♦ 5# 
in Air Density . 
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in a random variation in the^e parameters of about +10$. Air 
density, as described in equation 7> vas also given random variation# 
The multiplying coefficient, 1.22, of that equation vas allowed to 
vary normally with mean 1,2 2 and standard deviation of OCUU . This 
resulted in a random variation of air density of about +5/^ for any 
particular altitude. These variations seemed only to propagate into 
the final errors in x, Z, and V and into the extreme control 
excursions in angle of attack. Of these, only the variation in the 
final altitude causes concern. It varies between . h and 2.4 meters. 
Since the sink rate is about 3-24 metex s/second, it follows that 
a is corresponds to an error of about .6 seconds . time. The 
control system appears to handle these normal variations and 


uncertainties q* ite well. 



COHCLJSIOHS 


This paper describes cm application of linear optimal 
regulator theory to a nonlinear sinulation of an aircraft performing 
a helical landing approach. The nonlinear equations of motion are 
developed and are linearized (time varying coefficients) about a 
quasi-steady helical flight path. The nominal state time histories 
are given as explicit functions of time and a numerical method for 
determining the constant control inputs is presented. Control of the 
system to the nominal state trajectories is posed as a regulator problem 
with time varying weighting luatrices in the cost functional. A method 
of solving for the feedback gain matrix is reviewed. This theory is 
then implemented in a simulation of Boeing 73 7 -100, and system perfor- 
mance was measured for seven distinct approaches including flight 
path angles (desc nt angles) of three and six degrees. On each 
approacn , the aircraft was subject to large errors in initial values 
of state variables and to strong steady crosswinds. The system was 
also tested for sensitivity to normal variations in atmosphere and 
to reasonable uncertainties in the parametric description of the 
aircraft . Statistical data on 200 simulated landings is presented. 

The control system performed very well on all the three degree 
approaches and was reasonably insensitive to changes in the 
atmosphere and to parametric change? in the aircraft mo- 1. 

Performance, which was measured i ^ of terminal errors, 

violation of design constraints, and passenger comfort, was not 
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nearly as good on the six degree approaches. These approaches are 
more difficult to execute since descent rate and roll angle are 
'it led and the radius of the helix is halved. On the most 
difficult of these six degree approaches, the total flight time vas 
about sixty seconds, and performance vas definitely unacceptable. 
Evidence is presented to support the contention that total flight 
time is the critical factor and that performance on the six degree 
paths could be greatly improved by starting the approach from a 
higher altitude. This hypothesis needs to be further tested if 
six degree approaches are a requirement. 

The control system design method used in this study is relatively 
straight forward and is easy to implement with the aid of a modern 
computer. The only difficulty is in choosing the weighting matrices 
for the cost functional. It vas concluded in this paper that the 
weights should be time varying for the particular control task 
studied. This method should oe extended to the more complex simula- 
tion which includes both actuators and sensors and has six degree 
of freedom dynamics before being implemented on the aircraft. 
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APPENDIX A 


JACOBIAN MATRICES USED IN THE 
LINEARIZATION OF THE AIRCRAFT DYNAMICS 

The aircraft dynamics were linearized about a nominal 

state trajectory and nominal control time history. This was 

done in order to apply linear optimal regulator theory to the 

computation of the feedback gains. As is evident in equation 

(3) of the text, the Jacobian matrices 

A(t) * and B(t) = 

dx 3u 

are needed. They are 

0 0 0 

0 0 0 


3D 

da 


/m + (t/m) sin a 


3L 

T cos a + -gjj cos $ 


/mV (L/mV) sin $ 


- sin a/V 


sin <fr/(mV cos y) 


(L/mV) cos $/cos y 


" he " !l * 8 s \ ife ’ !? * s ^ "(no) (“ * °o> 


and 
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0 


cos Y cos 4» 



0 

0 

0 

0 

0 

0 


0 

0 

0 

0 

0 

0 


cos y sin 4» 


sin y 
2D/(mV) 

(T/m) sin a - g cos Y - (L/m) cos # ^ V 
(L/m) sin $ / (V 2 cos y) 


-2 


- V sin y cos 4* 

- V sin Y sin (p 

V cos Y 
g cos Y 
(g/V) sin Y 

(L/m) sin $ sin y / {V cos^ 


-V cos y sin 4» 
V cos y cos 4> 
0 
0 
0 

y) 0 


(A 2) 


An assumption that atmospheric density is constant over the range 
of altitude considered is incorporated into these equations. This 
assumption is discussed in the text. It is not made in the non- 
linear simulation which uses the feedback gains. 
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